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Abstract. The Darbroux transformation is generalized for time-dependent Hamilto- 
nian systems which include a term linear in momentum and a time-dependent mass. 
The formalism for the AT- fold application of the transformation is also established, and 
these formalisms are applied for a general quadratic system (a generalized harmonic 
oscillator) and a quadratic system with an inverse-square interaction up to A = 2. 
Among the new features found, it is shown, for the general quadratic system, that the 
shape of potential difference between the original system and the transformed system 
could oscillate according to a classical solution, which is related to the existence of 
coherent states in the system. 

PACS numbers: 03.65. Ge 03.65. Ca 
1. Introduction 

There has been great interest in using the Darboux transformation PU |2j for the analysis 
of physical systems and for finding new solvable systems. It has been shown that the 
transformation method is useful in finding soliton solutions of the integrable systems jHj 
and constructing supersymmetric quantum mechanical systems jlj. An excellent survey 
for developments and some applications of the transformation method is given in j^J. 

If a Darboux transformation is applied to a time-independent Schrodinger equation 
of confining potential, one of the immediate consequences is that the energy eigenvalues 
of the transformed system will be almost identical with those of the original system 
except for a finite number of addition(s) and/or deletion(s) to the spectrum 0. On the 
other hand, Abraham and Moses (AM) have developed an integral equation algorithm 
[Zj, which can also be used to find a new solvable system based on a known one. When 
applied to a harmonic oscillator, the AM algorithm gives a transformed system whose 
energy spectrum coincides with that of the harmonic oscillator except that the lowest 
eigenvalue has been removed. It has been shown that the result given by AM can also 
be derived through the factorization method [Hj and the factorization method has been 
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applied for various systems ||] . After all, the AM method and the factorization method 
are intimately related to the Darboux transformation [BJ E3 E] • 

In order to obtain a new solvable system by implementing the Darboux 
transformation, it is necessary to choose the auxiliary function of the transformation 
judiciously [0]. After the transformation being extended to include a time-dependent 
potential, the application of the Darboux transformation and the two-fold application 
of the transformation have been explicitly carried out for the harmonic oscillator model 
(with time-dependent frequency), to obtain new solvable systems |11[IT2]. For the simple 
harmonic oscillator, unphysical negative-energy eigenstates, which can be obtained by 
invoking the symmetry of the Hamiltonian, have been used as the auxiliary functions 
of the Darboux transformation. For the two-fold application of the transformation, a 
theorem has already been established in ^HJ for the choice of the auxiliary functions. 

In this article, we will generalize the Darbroux transformation to be applicable to 
time-dependent Hamiltonian systems which include a term linear in momentum and a 
time-dependent mass. The formalism for the A-fold application of the transformation 
will also be established. The transformation method will then be applied for a general 
quadratic system (a generalized harmonic oscillator) and a quadratic system with an 
inverse- square interaction up to N = 2. The general quadratic system is known to 
be related to a simple harmonic oscillator through unitary transformations, which may 
be a reason for the existence of coherent states in the simple harmonic oscillator [15J. 
Even for the cases already considered [TT1 IT2~] , new features of the Darboux transformed 
systems will be given. In particular, it will be shown that, the shape of the potential 
difference between the original Hamiltonian and the transformed one could oscillate 
according to the classical solution of the quadratic system, which is related to the 
existence of coherent states in the system. 

2. A generalization of the Darboux transformation 

For smooth functions P(x), Q(x), we assume a function u(x) satisfies 

^ + P(^ + (Q(x)+C)u = 0, (1) 
where C is a constant. If <p(x) satisfies the linear differential equation 

' / %P(*)^ + Q(*)0 = O, (2) 



dx doc 
Darboux shows that the following equation is true [T] 

d 2 „ . , d dP „ , s , d 2 In u . 



dx 2 ^ ^ dx dx ^ ^ v dx 2 




(3) 



This Darboux transformation has immediate consequences for a time-independent 
Schrodinger equation jHj. 
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For the extension of the transformation to be applicable for the time-dependent 
system described by the Hamiltonian 

2 

H = JLr-z + (R(x, t)p + P R(x, t)) + V(x, t), (4) 



2M(t) 

we consider the operator 

0(t,x)= - iti^ + H 



t, 2 



m mt) §?-*M*,t)± + K( x ,t)} + v( x ,t). (5) 

In H, M(t) denotes time-dependent mass and the term proportional to R(x, t) is 
included to resemble a three dimensional model which is under a vector potential. A 
wave function of the system of H should satisfy the time-dependent Schrodinger equation 

0(t,x)ijj(t,x) = 0. (6) 

We also introduce the auxiliary functions Uk satisfying 

0(t,x)u k {t,x) = (k = l,2,---,N). (7) 

We assume that Uk is a smooth function for any finite x, but Uk need not be square- 
integrable as it is just an auxiliary function. 

From (6) and (7), it is straightforward to show that 

h 2 d 2 \n Ul 



0(t,x) -2ihR'(x,t) 



4> ld (x,t) = o, 



M{t) dx 2 



where 



^ ld (x,t)=U'-^\, 



(8) 



(9) 



and ' denotes the partial derivative with respect to x. Equation (8) is a generalization 
of the Darboux transformation. 

The Darboux transformation can be applied repeatedly. For the description of 
fc-fold transformations, we define the Wronskian determinants as 



Ui 

u' 



U 2 



(fc-1) (fe-1) 



u 



(fe-1) 



I 

1 



u 



u 2 

I 

2 



u 



(fe) (fe) 



u k 

I 

k 



u 



u 



(fe) 



(10) 
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Making use of Crum's formula [2] 

= W k —W k ^ - W k ^—W k , (11) 



one can find that ip (x,t) satisfying 



h 2 d 2 \nW k \ M 



0(t, x) - 2ikhR'(x, t) - ^ J ^ d (x, t) = (12) 



is given by 



V> W (*,t) = %* (13) 
yVk 

There is a difficulty in interpreting equation (12) as a Shrodinger equation, since 
the associated Hamiltonian is, in general, not Hermitian. As discussed in special cases 
[TT] . this difficulty can be resolved if there exists a purely time-dependent function a k (t) 
satisfying 

-\na k = 2kR> - i-^^ In (14) 
dt 2M(t)dx 2 W k y ! 

where W k denotes the complex conjugate of W k . In this case, equation (12) can be 
rewritten as 

where 

i// eD (x,t) = a k {t)i/> kd (x,t). (16) 

If V(x,t) is a real function, equation (15) is the Schrodinger equation of a system 
described by the Hermitian Hamiltonian 

a -= 2m + <*<*• t)p + ()) + v(x ' t] ~ ^m^r- < 17 > 

One of the crucial conditions for ip kD (x, t) being square-integrable is that W k should not 
have any zero in the entire space of x. Indeed, the derivations of (12) and (13) are not 
valid for the zeros of Wi, W 2) • • • , W k , while, if W k has no zero, the formulas would still 
be useful even for the cases that Wi, W^, ■ ■ ■ , W k -\ have zeros, as examples will show 
later. 

In addition to ip hD ( x ,t), other solutions of the Schrodinger equation of H k would 
exist For the H\ system, from the fact that 

8 1 H 2 ( d 2 1\ +„(dl\ / _ _ _1 lh 2 (d 2 \nu\ , . 

a solution ip]^ satisfying [—ih(d/dt) + H^ipl® = is given as P ITT] 

< D M = — = 4r- ( 19 ) 
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3. A general quadratic system 

For an application of the transformation of the previous section, it may be essential to 
find U\ or Wk which does not vanish over the whole coordinate space. Moreover, for an 
Hermitian Hamiltonian, such auxiliary solutions are required to support the existence 
of a purely time- dependent function, defined in (14). For the non- vanishing 

U\ of the harmonic oscillator, unphysical negative-energy eigenstates can be used as 
the auxiliary functions for the transformation. For the two-fold application, the sign 
theorem established in may be used to find non- vanishing Wi. Nevertheless, it has 
been shown that the application and two-fold application of the Darboux transformation 
are possible for the harmonic oscillator (of time-dependent frquency) with and without 
an inverse-square potential [TTJ [T2] . 

In this section, we will show that the application and two-fold application of the 
Darboux transformation are possible for a general quadratic system. Even for the cases 
already considered [HI [T2] , new features of the applications will be found. 



3.1. Darboux transformation 

A general quadratic system can be conveniently described by the Lagrangian [1] 
L Q = \m{1)x 2 - l -M(t)w 2 (t)x 2 + F{t)x + j t (M{t)a(t)x 2 ) + ^(6(t)x) + f(t), (20) 

where the overdot denotes the derivative with respect to t. This Lagrangian gives the 
classical equation of motion 

4- (M(t)x) + M(t)w 2 (t)x = F(t) (21) 
at 

which is that of a generalized harmonic oscillator of mass M(t) and frequency w(t) in 
an external force F(t). The corresponding Hamiltonian is written as 

H Q = - a(t)[px + xp] + l -M{t)c{t)x 2 - M.p + d(t)x + (^ - /(*)), (22) 

where 

c(t) = w 2 + 4a 2 - 2d - 2^a, d(t) = 2ab - b - F. (23) 

We assume that Hq is Hermitian. The general solution of equation (21) is a linear 
combination of a particular solution x p (t) and two linearly independent homogeneous 
solutions u(t),v(t). We assume x p (t),u(t),v(t) are real, and define pit) and a time- 
constant Q, for later use, as 

p = ^u 2 [t) +v 2 (t), VI = M(t)[v(t)u(t) - u(t)v(t)\. (24) 

Since the general quadratic system can be obtained from a simple harmonic 
oscillator through unitary transformations, we first consider the simple harmonic 
oscillator of unit mass and frequency. The Schrodinger equation for the simple harmonic 
oscillator, whose time is r, is written as 

d h 2 d 2 x 2 
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Making use of the invariance of the Schrodinger equation under the exchange of r <-» — r 
and x <r-> ix, from the well-known wave functions 



-i(n+l/2)r 



exp 



X 



x 



(2 n nWTihYI 2 

one can find auxiliary functions satisfying the same Schrodinger equation as 



<(r, x) 



e i(n+l/2)r exp 



X 

2h 



Hr, 



IX 

7h y 



n 



0,1,2,.-.. 



(26) 



(27) 



Indeed, v^(r,x) for even n has no zero over the whole coordinate space, and it can be 
used as an auxiliary function of the Darboux transformation. In this subsection, we 
restrict our attention on the cases that n — 0, 2, 4, . . .. 
By defining 



d 

O s (t,x) = -ih— + H s , Q (t,x) 



•n 9 



Hq, 



and if r and t are related by 



dr 



-dt. 



M(t)p\t) 

then from the results given in ^2] one can find that 



Q (t,x) = U G U F U s O s (r,x)UlulU G \ T=T{t) 



where 



G = exp[-(M(t)a(t)x 2 + b(t)x + j f(z)dz)} 
U F (x p ) = exp[^(Mx p x + 5(x p ))}exp(-^x p p) 



u s ( P ,n) 

with 5 defined through the relation 

6{x p ) 



exp[^M^x 2 ] ex Ph^ ln ( fy( x P + P x )} 



\mw 2 x\ - X -Mx\. 



(28) 
(29) 

(30) 

(31) 
(32) 
(33) 

(34) 



A point that should be mentioned is the unitary operators Uf and Us are not 
unique since they depend on the choice of classical solutions. Instead of {u,v} and x p , 
one can take another set of two linearly independent homogeneous solutions {u, v} and 
a particular solution x p of equation (21). After defining p, Q and 5 from {u, v} and x p 
as p, Q and S are defined from {u, v} and x p , one can find that the unitary relation (30) 
is also valid with the unitary operators Uf(x p ), Us(p, fi). 

Making use of the unitary relation, the normalized wave functions of the system of 
Hq have been given in pH] as 



n 



V2 m m\ \7chJ ^(t) 



u(t) - iv(t) 



exp 



x exp 



l r 



M(t)a(t)x 2 + (M(t)x p (t) + b(t))x 
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x exp 
x H m 



(x - x p {t)f 

2h 

iQx — x p (t) 



+ lM(i) ^. 



a pW (35) 

The unitary relation can also be used to find a solution v®{t,x) of the Schrodinger 
equation 



as 



O(t,x)v%(t,x)=0, 
v%(t,x) = U G U F (x p )U s (p,n)v s n (r,x)\ T=T{t) 



(36) 



x exp 
x exp 



u(t) + iv(t) 

~~W) 



exp 



h 



6(t) + / f(z)dz 



M{t)a{t)x 2 + (M{t)x p {t) + &(*)) 



1 

h 



2ft 



i 



. JVtx — x p (t) 

" n p(t) 



• (37) 



\P 2 (t) w />(*). 

From the properties of the unitary transformations, it is manifest that v®(t,x) for even 
n does not have a zero over whole coordinate space. 

Since, for even n, H n {^\[^ x ~p(Jf' ^j is a real function, it is easy to find that 



2M or 2 p 



(38) 



which shows, for the choice of uf = = v®, that a®{t) defined in (14) can be found 



as 



(39) 



up to a normalization constant. Therefore, the transformation formalism developed in 
the previous section can be applied with uf(t,x) = v^(t,x) to give the solvable model 
described by the Hamiltonian 

2" 



Hf(t,x,p) = H Q - J ^ + 4n 



hn 



hn 



Mp 



n2 



(n-1) 



H n - 2 (z) 
H n {z) 



n 



. H n (z) 



(40) 



where z = iyj^-^. If we adopt the notation that H_ 2 (z) = H_ 1 {z) = 0, equation (40) 



is valid for n = 0, 2, 4, • • •. The magnitude of AV (= #r Q -H Q + j£k) is O(ft). Since 



AV^ vanishes in the limit of \z\ — > oo, AV/h approaches except for the region in the 
vicinity of x p where the width of the region is 0(\fU). 

From equations (10), (13) and (16), the unnormalized wave functions tp^(t,x) 
satisfying the Schrodinger equation 



(41) 
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can be given as 



C Q M= ,50 m [~M(s p - x p ) + 



+ 2mi 



r 

2in\ 



n 



'fUtf„-l(z)- 



(42) 



where u> 



ft p H m (w) ""V ft p if n (,z) 

o £-^p ^ g ^ e general quadratic system [HS , a different choice of 
the classical solutions u(t),v(t) x p (t) gives different wave functions. When we choose 
u(t) = u(t), v(t) = v(t), and x p (t) = x p (t), for even integer n, an unnormalized tp^(t, x) 
is written as 



vm + l 



m + if 

, P , 



r- .H n -i(iw) 
Vv+i + V2m- 



HJiw) 



(43) 



for m = 0, 1,2, For the case of M(t) = 1 and a(t) = 0, i[)^(t,x) in equation (43) 

reproduces the wave functions given in up to normalization. From equation (19), 
another wave function satisfying (— ik-j^ + ip™® = is given as 

tfS = "4- (44) 

pUn 

For non- negative even integer n, it is clear that i/iffi (m = —n, 0, 1, 2 . . .) are square- 
integrable. 

Even for the simple harmonic oscillator of unit mass and frequency, since the shape 
of the probability density of a wave function could breathe and oscillate [HJ [T5] , the 
region, where AV/h is different from by a certain amount, could breathe and oscillate. 
If we take the classical solutions as u(t) = u(t) = cost, v(t) = v(t) = csint (c ^ 0), and 
x p (t) = dcost, the breathing and/or oscillating behavior of the non- vanishing region 
appears when c / 1 and/or d ^ 0, respectively. For the choice c = 1 and d — 0, H™® 
becomes 



b p 2 +x 2 ) _ h _ h 2^ hlHn(i _ X 

2 ox* 



ny Wh> 



(45) 



Through the same choice of classical solutions, for non-negative even integer n, one can 
also find the eigenfunctions of H™ s as 



m- 



V n 



(46) 



whose eigenvalues are (m + |)ft and — (n + |)ft, respectively, as expected in p. 



3.1?. ^4 two-fold transformation 

In this subsection, we will show that the two-fold application of the Darboux 
transformation is also possible for a general quadratic system. It will also be manifested 
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that, in the two- fold transformation, the breathing and oscillating behavior would still 
appear in the transformed systems of a simple harmonic oscillator. 
By taking 

ui(t,a;) = i> n {t,x), u 2 (t,x) = ip n+1 (t,x), (47) 
one can find that 

. o m i 1 r 

^^ 2 {t,x)J n {w), (48) 
where 

d d 
J n {w) = H n (w)—H n+1 (w) - H n+1 {w)—H n {w) (49) 
aw dw 

with w = sf^^-- Making use of a recursive relation among the Hermite polynomials 
[TT)] . one can show that J n {w) = 2H 2 (w) + 2nJ n -\(w). This relation, with the fact 
Jo = 2, shows that J n {w) is positive definite for all x, which in fact has been implied in 
[15] . One can also find that 

-i4!U^ = 4a + 2^. (50) 
2Mdx 2 W 2 p K ' 

Equation (50) shows that a 2 can be given in accordance with equation (14) as 

a Q 2 (t)= P 2 (t) } (51) 

up to a multiplicative constant. 

With the W 2 , the two-fold transformation can, therefore, be applied to the quadratic 
system to give the transformed Hamiltonian 

^= H ° +2 w-^ Mw) - (52) 

For the case of M(t) = 1 and a(t) = 0, reduces to the Hamiltonian found in |T2]. 
Making use of equations (10), (13) and (16), it is also possible to find the wave functions 
of the system of H^- F° r non- vanishing W 2 , it turns out that the same set of classical 
solutions must be used for both u\ and u 2 in equation (47). However, in obtaining wave 
functions, a different set of classical solution could be used for ip in equation (10), as in 
the Darboux transformation of the model. 

4. A quadratic system with an inverse-square interaction 

In this section, we will consider the application of the Darboux transformation for the 
quadratic system with an inverse-square interaction described by the Hamiltonian 

Hm = 2M{t) ~ a{t){xP +PX) + \ M[t ^ x2 + W(t)^ ^ 
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denned on the half line x > 0, where g is a constant. The system of H in is related to 
the system described by the Hamitonian 

W = V - + -x 2 + °- 

through a unitary relation [15] . If u(t),v(t) denote the homogeneous solutions of 
equation (21) as in the previous section, and pit) and Q are defined by equation (24), a 
wave function of the system of Hi n is given, for non-negative integer n, as 




x exp 



(2n+a+l) 



fix 2 



'fix 

hp 2 



2 \ (2a+l)/4 



hp 2 



(54) 



,P P 

where a is defined through the relation g = \(oi+ |)(a — |)^ 2 , and L° is the Laguerre 
polynomial For a > —1, i[)™(t,x) is square-integrable on the half line. 

Through similar procedures used in the previous section, one can find the auxiliary 
function 

(2n+o+l) /n 2 x (2o+l)/4 



1 

Vp V 

x exp 



/ u + if ' 



x 

2h 



satisfying 




^ + iM^ + 2iMa 

,P P 



0. 



Q x 2 
hp 2 



(55) 



(56) 



From the properties of the Laguerre polynomial, v™ has no zero on the half line, and it 
is easy to see that 

. h d 2 . vl r . 1 n p 
-lnd?- = 2a + -. 

P 



(57) 



2M dx 2 v™ 

Equations (14) and (57) show that a™{t) can be given as 

<(t) = p{t) (58) 

up to a normalization constant. The Darboux transformation can therefore be carried 
out, with u\ = v™(t, x), to find the Hamiltonian of a solvable model as 

h 2 <9 2 (<x n ) 



(59) 



2M dx 2 

In implementing the formulas (10), (13) and (16), different sets of homogeneous 
solutions of equation (21) can be used for v™{= U\) and V'm( = to find the general 
expression of the wave functions of the JJ™ ,m system, as in the general quadratic system 
of the previous section. For simplicity, however, we only consider the case that the same 
set {u(t),v(t)} of the homogeneous solutions is used in t>™ and ip™. The wave functions 
•0™'*" of the Hi' m system, if we adopt the notation L_i(y) = 0, is given as 

2 
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x 



(m — n— l)y + (n + a 



-y) 



K(-y) 



in 
m ' 



for m — 1, 2, 3, . . ., where 
Qx 2 

y 



hp 2 



For m = 0, the wave function of the ' system is given as 



2^ 



-(n + l)y + (n + a) 



-y) 



fo- 



il 

(60) 
(61) 
(62) 



L%(-y) 

It is clear that ijj^ m (t, x) is square- integrable on the half line for a > 0. Though another 
formal solution of the Schrodinger equation of the H™ ,m system can be found through 
equation (19), that solution turns out to be not square-integrable for a > 0. 

As implied in ^3] , for the two-fold application of the Darboux transformation, the 
auxiliary functions can be chosen as 



ui{t,x) = fon(t,x), u 2 (t,x) = Cn(M), 



which gives 



where 



111 



n + 1 



u — w 



4n+2 



(n + a)\ V n + 1 



a: 



2 Qx k ( \ 

j- 2 K n (y), 



K n (y) = L«{y)^-L a n+1 {y) - L a n+1 (y)^(y). 



(63) 



(64) 



(65) 



Making use of a recurrence relation among the Laguerre polynomials, one can easily 
find that 



K n (y) 



1 



(66) 



n + 1 v " w " ' n + 1 
With the fact that K = —1, (66) proves that K n (y) is negative-definite for all n, so 



that W™ does not vanish for all < x < oo. The fact 



h d 2 



— i 



In 



2M dx 2 Wl n 



4a + 2 



P 

J 

P 



(67) 



shows that the Hermiticity condition (14) of the Hamiltonian can be satisfied with 

a?(t)=p 2 (t). (6 
With WjfN from (17), the Hamiltonian of the transformed system can be found as 

h 2 d 2 



2h 2 (a + i) nn 

tin n H ——^ r 



\nK n (y). 



(69) 



Mx 2 Mp 2 M dx 2 

It may be possible, through the formulas (10), (13) and (16), to find the wave functions 
of the system of H™. 
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5. Discussion 

We have generalized the Darboux transformation to be applicable to a general 
one-dimensional time-dependent Hamiltonian system. The formalism for an iV-fold 
application of the transformation has also been established. It has been shown that an 
Hermitian system can be found from the transformed system, if the auxiliary function(s) 
of the transformation could support the existence of a purely time-dependent function 
satisfying a certain condition (equation (14)). The formalisms have been applied to a 
general quadratic system and a quadratic system with an inverse-square interaction. As 
the potential difference between the original and the transformed systems is calculated 
from a (formal) solution of the Schrodinger equation, due to unitary relations responsible 
for the existence of coherent states in the system, the shape of the potential difference 
could oscillate according to a classical solution for the general quadratic system. 

Since the potential difference is calculated from a solution of the Schrodinger 
equation, the difference depends on the Planck constant Ti. In the cases considered, 
the magnitude of the potential difference is 0(h). The range within which the potential 
difference effectively depends on the space coordinate would also be described in terms 
of h, as has been shown explicitly in the examples considered. These features are not 
usual in the standard text book of quantum mechanics, and the potential difference 
appears as if it is a quantum correction which depends on the position. 

It should also be mentioned that, for the systems considered in this paper, 
there may exist other choices of auxiliary functions which lead to Hermitian systems. 
Only for limited cases, however, it may be possible to find exact solutions of the 
Schrodinger equation of a time-dependent system. Even if the solutions are found, the 
auxiliary function should be chosen judiciously so that the transformed system could 
be Hermitian. As an example, vf + v ® is a solution of the Schrodinger equation of the 
general quadratic system, so that the Darboux transformation can be formally applied 
to it. However, with this solution, one can not find a{t) satisfying (14). 
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